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Abstract. An exact analytic solution is obtained for a uniformly expanding, neutral, highly conducting 
plasma sphere in an ambient dipole magnetic field with an arbitrary orientation of the dipole moment in 
the space. Based on this solution the electrodynamical aspects related to the emission and transformation 
of energy have been considered. In order to highlight the effect of the orientation of the dipole moment in 
the space we compare our results obtained for parallel orientation with those for transversal orientation. 
The results obtained can be used to treat qualitatively experimental and simulation data, and several 
phenomena of astrophysical and laboratory significance. 

PACS. 03.50.De Classical electromagnetism, Maxwell equations - 41.20.Gz Magnetostatics; magnetic 
shielding, magnetic induction, boundary-value problems - 41.20.Jb Electromagnetic wave propagation, 
radiowave propagation - 52.30.-q Plasma dynamics and flow 



1 Introduction 

There are many processes in physics which require the 
solution of boundary and initial value problems. When 
the boundary is immobile the standard techniques for the 
solution of the problem are well known. However, the 
situation becomes very complicated introducing a mov- 
ing boundary into the problem. Usually this excludes the 
achievement of an exact solution of the problem and the 
development of some approximate methods is desirable 
[IH3]. The simplest case often allowing an analytical so- 
lution is the moving plane boundary. With the help of a 
time-dependent transformation it is possible to immobilize 
the boundary but increasing essentially the complexity of 
the underlying dynamical equations. The purpose of this 
work is to present an example of an exactly solvable mov- 
ing boundary and initial value problem and generalize the 
results of our previous paper [4] . 

The problems with moving boundary are very impor- 
tant in fundamental mathematical physics but they have 
also important applications in many area of physics. For 
instance, an important example is the expansion of a plasma 
cloud with a sharp boundary across an ambient magnetic 
field. Such kind of processes are particularly of interest 
for many space and laboratory investigations (see 5 9 
and references therein). This is a topic of intense interest 
across a wide variety of disciplines within plasma physics, 
with applications to laser generated plasmas [5] , solar [TU] 
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and magnetospheric [TT1IT2] physics, astrophysics [13] , and 
pellet injection for tokamak refueling |14j . 

The dynamics of the plasma cloud has been studied 
both numerically [5HE] and analytically [T51 - [T9"] . Usually 
the plasma is treated as a highly conducting media shielded 
from the penetration of the magnetic field. Thus the mag- 
netic field is zero inside. An exact analytic solution for a 
uniformly expanding plasma sphere in an external uni- 
form magnetic field has been obtained in [T51Q35] both 
for non-relativistic |18] and relativistic |19] expansions. 
Within one-dimensional geometry a similar problem has 
been considered in [T5] . In the previous papers [3] and [T7] 
we obtained the exact solutions for the uniform relativistic 
expansion of the plasma sphere and cylinder in the pres- 
ence of a dipole and homogeneous magnetic fields, respec- 
tively. However, in Ref. [3] we have considered a somewhat 
simplified situation when the dipole moment is directed 
to the centre of the plasma sphere thus providing an az- 
imuthal symmetry to the problem. In the present paper we 
generalize the solution obtained in [5] considering a simi- 
lar problem of the expansion of the plasma sphere in the 
presence of a dipole magnetic field but with an arbitrary 
orientation of the dipole moment in the space. 



2 Moving boundary and initial value problem 

We consider the moving boundary and initial value prob- 
lem of the relativistic (v < c, where v is the radial velocity 
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of the sphere) expansion of the plasma sphere in the vac- 
uum in the presence of a dipole magnetic field. Consider 
the dipole with the magnetic moment p and a neutral 
infinitely conducting plasma sphere with radius R(t) lo- 
cated at the origin of the coordinate system. The dipole 
is placed in the position ro from the centre of the sphere 
(R(t) < ro). The orientation of the dipole moment is given 
by the angle i? between the vectors p and ro. In contrast 
to our previous paper [1] we consider here an arbitrary 
orientation of the dipole moment in space with arbitrary 
The plasma sphere has expanded at t = to its present 
state from a point source located at the point r = 0. It 
is convenient to introduce the vector potential of the in- 
duced and the dipole magnetic fields. The magnetic field 
of the dipole is given by Ho = V x Ao, where the vector 
potential A is 



A = 



p x (r - r ) 



(1) 



As the spherical plasma cloud expands it both per- 
turbs the external magnetic field and generates an electric 
field. Within the spherical plasma region there is neither 
an electric field nor a magnetic field. We shall obtain an 
analytic solution of the electromagnetic field configuration 
for an arbitrary 

When an infinitely conducting plasma sphere is intro- 
duced into a background magnetic field the plasma ex- 
pands and excludes the background magnetic field to form 
a magnetic cavity. The magnetic energy of the dipole in 
the excluded volume, i.e., in the volume of the conducting 
sphere is calculated as 



8tt 



dr = Q(£) = Q» (0 cos 2 + Q ± (f ) sin 2 0, (2) 



derived in Appendix A of Ref. [I] one obtaines: 
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where Pjf (x) is the generalized Legendre polynomials with 
v = 0, 1 and P e (x) = P${x). 

Since the external region of the plasma sphere is devoid 
of free charge density, a suitable gauge V-A = allows the 
electric and magnetic fields to be derived from the vector 
potential A. Having in mind the symmetry of the problem 
with respect to the a;z-plane it is sufficient to choose the 
vector potential in the form A r — An r , and 

oo 

A B = A oe + sin sin <p At (r, t ) P[ (cos 6) , (6) 
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and the components of the electromagnetic field are given 
by E r = and 



where £ = R/rn < 1, Qr is the volume of the plasma 
sphere, Q,|(0 = (j? / '2rg)£/i(0 > Q±(0 = (P 2 /^f 2 (0 
and the functions fi(x) and faix) are determined in Ap- 
pendix A by Eqs. (A. 3) and (A. 4), respectively. Here Qu (^) 
and Q±(£) are the magnetic energies escaped from the 
plasma volume at parallel (i9 = 0) and transversal (i? = 
7r/2) orientations of the dipole, respectively. The energy 
Q(£) increases with decreasing and reach its maximum 
value at $ = or $ = tt, that is the magnetic moment p 
is parallel or antiparallel to the symmetry axis ro. 

We now turn to solve the boundary problem and cal- 
culate the induced electromagnetic fields which arise near 
the surface of the plasma sphere due to the dipole mag- 
netic field. Since the sphere is highly-conducting the elec- 
tromagnetic fields vanish inside the sphere. In addition, 
the normal component of the magnetic field H r vanishes 
on the surface of the sphere. To solve the boundary prob- 
lem we introduce the spherical coordinate system r, 9, tp 
with the z-axis along the vector ro and the azimuthal an- 
gle ip is counted from the plane (irz-plane) containing the 
vectors ro and p. The spherical coordinate 9 is the angle 
between the radius vector r and ro. The vector potential 
([3} at r < ro can alternatively be represented by the sum 
of Legendre polynomials. Using the summation formulas 
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r dr 
IdAe 



dt 



(8) 

(9) 
(10) 
(11) 



In Eqs. ([6]) and (0 the prime indicates the derivative with 
respect to the argument and A.e(r,t), Bi(r,t) and Ce(r,t) 
are the expansion coefficients which are determined from 
the boundary and initial conditions. It is straightforward 
to show that the gauge V • A = is satisfied automatically 
if Ci(r,t) — Ai(r,t). Previously in Ref. [1] we have con- 
sidered the case of the parallel or antiparallel orientation 
of the dipole when i? = or = tt. In this case Eqs. ©- 
(ITTj) are essentially simplified because A r = A 0r = 0, 
Ag = A oe = and H v = E r = Eg = 0. Then the problem 
is reduced to the evaluation of the vector potential A v 
which involves only the coefficient Bg(r,t). 

For arbitrary ■& the equations for the expansion coeffi- 
cients At(r, t) and Bi(r, t) are obtained from the Maxwell's 
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equations 



A(r, t) — A (r) at r ^ ct and A(r, t) — at r < vi. Here 



d 2 A e + 2 d At 
dr 2 r dr 



£(£ + !) 
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1 d 2 A e 
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(12) pe ( z ) = 2 i £l(z 2 



Similar equation is obtained for the quantity Be (r, t) . This 
equation is to be solved in the external region r > R(t) 
subject to the boundary and initial conditions. Here R(t) 
is the plasma sphere radius at the time t. The initial con- 
ditions at t = arc 



Ai(r,0) = B e {r,0) = 0, 



dAi(r,0) _ dB e (r,0) 



dt 



dt 



= 0. 



(13) 

The first initial condition states that the initial value of 
A is that of a dipole magnetic field, A(r, 0) = Ao(r). 
The second initial condition states that there is no initial 
electric field. Boundary conditions should be imposed at 
the spherical surface r = R(t) and at infinity. Because of 
the finite propagation velocity of the perturbed electro- 
magnetic field the magnetic field at infinity will remain 
undisturbed for all finite times. Thus, for all finite times 
Ae(r,t) — s> and Be(r,t) -> at r -> oo. The boundary 
condition at the expanding spherical surface is H r — or, 
alternatively, 



A t (R,t) 



' f - \ , B t (R,t) 



£+1 \r 



We consider the case of the uniform expansion of the 
plasma sphere R(t) — vt with a constant expansion ve- 
locity v. This special case of the uniform expansion falls 
within the conical flow techniques which has been applied 
previously in Refs. |17H19j . From symmetry considerations 
one seeks a solution for the total (i.e., the unperturbed 
potential Ao(r) plus the induced one) vector potential of 
the form Ai{r,t) = r v W t {Q and B e (r,t) = r v $ t {C) with 
( = r/ct, where c is the velocity of light. Here &e.(Q and 
<Pe(C) are some unknown functions and v > 0. Having in 
mind the symmetry of the unperturbed magnetic field and 
also the boundary conditions (IT4l) it is sufficient to choose 
the parameter v as v = £. The equation for the vector po- 
tential A(r,i) is obtained from the Maxwell's equations 
which for the unknown functions ^eiC) an d &e(C) yields 
the same ordinary differential equation 

C(i - C 2 )^'(C) + 2(* + 1 - C 2 )^(C) = o. (15) 

This equation is to be solved in the external region r > 
R(t) subject to the boundary and initial conditions. The 
boundary condition at the expanding spherical surface 
is H r — which is equivalent to the relations \&i(0) = 

-m& and = ~& ( see Ec i' ®) with P = 

v/c < 1. In addition, imposing that A(r,i) = Ao(r) at 
r ^ ct we obtain another boundary conditions ^(1) = 
<£ £ (1) = 0. Thus, the solution of Eq. (JTS]) subject to the 
initial and boundary conditions may be finally written in 
the form at vt < r < ct 



A e (r,t) = —B e (r,t), B t {r,t) 



rYptil/O 



ro 



(16) 



1) + P," 



\Z)= / (T 2 ~lYdT (17) 



and P£(z) are the generalized Legendre functions with 
z > 1, [i — £, and v = — I — 1. 

The electromagnetic field components are determined 
according to Eqs. ©-(ll])- From Eqs. ©-(HU and (Tl6l) 
it is straightforward to show that the boundary condi- 
tions on the moving surface (see, e.g., [31[2D] for details), 
E(R) = -i[vxH(i?)] (or E e (R) = PH V (R) and E V (R) = 
—f3Hg(R)), are satisfied automatically. 

It should be emphasized that all above results are valid 
only for R(t) < ro or t < r^/v. At the time t = r$/v the 
plasma sphere reaches to the dipole which will be com- 
pletely shielded by the plasma. Therefore at t ^ r^/v the 
total electromagnetic field vanishes. 

Consider now briefly the non-relativistic limit of Eqs. © 
and ([7]) with the expansion coefficients (|T6"| recalling that 
Ci(r, t) = Ag(r, t). This limit can be obtained using at £ — » 
and j3 — > the asymptotic expression pi{l/C,)/pt{l/ (3) = 
(£/C) 2£+1 = {R/r) 2l+1 (see, e.g., Ref. [21] for the asymp- 
totic expansion of the generalized Legendre functions P^{z) 
with z > 1) which yields 



A g (r, t) = A oe (r) sin -d sin tp 

r R 

OO 

^(M)=^(r)-^E(7) 



(18) 



£+1 



(19) 



1 d 

cos -dP} (cos 8) + sin cos <x>- tttP} (cos ( 

ey ' £ + 1 dd e 



where r* = i? 2 /ro < R. The remaining ^-summations 
in these expressions can be performed using the sum- 
mation formulas derived in Ref. [1]. As has been shown 
in [4] the resulting induced vector potential in the non- 
relativistic limit represents a sum of the dipole and the 
quadrupole terms. These fictitious dipole and quadrupole 
('image' dipole or quadrupole) are located in the rrz-plane 
inside the plasma sphere at the distance r» = (R 2 /r 2 )ro 
from the centre. Finally, in the lowest order with respect 
to the factor /? the components of the electric field are 
given by Eg = —/3dAg/dR, E v = —fidA^/dR, where Ag 
and A v are determined by Eqs. (fT5)) and (IT^I . respectively. 



3 Energy balance 

In the problem of the plasma expansion in an ambient 
magnetic field it is important to study the fraction of en- 
ergy emitted and lost in the form of electromagnetic pulse 
propagating outward of the expanding plasma [41 IT5UT8] . 
In this section we consider the energy balance during the 
plasma sphere expansion in the presence of the magnetic 
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dipole with arbitrary orientation thus generalizing our pre- 
vious treatment [4] with specific orientations ■& — 0, it. It is 
assumed that the plasma sphere of the zero initial radius 
is created at t = 0. When it starts expanding, ambient 
dipole magnetic field Ho is perturbed by the electromag- 
netic pulse, H' = H — Ho, E, propagating outward with 
the speed of light. The tail of this pulse coincides with the 
moving plasma boundary r = R(t) while the leading edge 
is at the information sphere, r = ct. Ahead of the leading 
edge, the magnetic held is not perturbed and equals Ho 
while the electric field is absent. 

For the evaluation of the energy balance we employ 
the Poynting equation 



pressure and the decrease of the electromagnetic energy 
in a control volume f2 c . 

Consider now explicitly each energy component Ws(t), 
Wj(t) and AW E M(t)- W s {t) is evaluated from Eq. (gig) . 
In the first expression of Eq. (|2"2"j) the ^'-integral must be 
performed at *f ^ t' < t (t < £*) since at < t' < ?f the 
electromagnetic pulse does not reach to the control surface 
yet and S r (r c ) = 0. Using the summation formulas derived 
in Appendix from Eqs. ([5])-(fIT]). (fT6"|) and we obtain 



VS 



-j-E 



d E 2 + H 2 
dt 8^ 
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where S 



[E 



Jee g 
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j v e v (with 



x H 



is the Poynting vector and j = wnere r)-r c /ct<l, and 



= leu, I = 1) is the surface current 



density. The energy components involved in Eq. (|20p have 
been evaluated in detail in [3] . We recall here some results 
for completeness. In order to calculate the energy emitted 
to infinity we should integrate the Poynting vector over 
time and over the surface S c of the sphere with radius 
r c < (control sphere) and the volume fl c enclosing the 
plasma sphere (r c > R or ^ t < r c /v). Integrating over 
time and over the volume f2 c Eq. (|20[) is represented as 



where 



W s (t) = Wj(t) + AW EM (t), 

pt p 7T />2tT 

W s {t) = r 2 c / dt' / smOdO / S r d<p, 
Jo Jo Jo 

Wj{t) = - f dt' f j.Edr. 

Jo Jn a 



(21) 



(22) 



Here S r = ^(EgH^ 



is the radial component 



of the Poynting vector. Also Wem(*) an d ^Wem(£) 
W / em(0) — Wem(*) are the total electromagnetic energy 
and its change (with minus sign) in a volume fi c , respec- 
tively. Wj(t) is the energy transferred from plasma sphere 
to electromagnetic field and is the mechanical work with 
minus sign performed by the plasma on the external elec- 
tromagnetic pressure. At t = the electromagnetic fields 
are given by H(r, 0) = H (r) and E(r, 0) = 0. Hence 
Wem(0) is the energy of the dipole magnetic field in a 
volume fi c and can be calculated from Eq. @ by replac- 
ing R by r c . Thus Wem(0) = Q(u), where u — r c /ro < 1. 
Taking into account that H = E = in a plasma sphere, 
the change of the electromagnetic energy Z\Wem(0 m a 
volume f2 c can be represented in the form [4] 



AW EM (t) = Q(u) 



E 2 + H 2 



Szr 



(23) 



where fl' c is the volume of the control sphere excluding 
the volume of the plasma sphere. Hence the total energy 
flux, Ws(t) given by Eq. (f2"2"j) is calculated as a sum of the 
energy loss by plasma due to the external electromagnetic 



a*(i?) = (£+!) cos 2 d + - sin 2 i? 



,(24) 



(25) 



is the orientation factor of the dipole magnetic field. In 
the non-relativistic limit, /3 — > 0, using the asymptotic 
expression (see, e.g., Ref. 121 ) pt(z) — z 2e+1 /(2£+ 1) at 
z — > oo, as well as the summation formulas of Appendix 
from Eq. (j2~4"]) we obtain 



W s (t) 

„2 



+ 



P 



2Q(0 

^•3 



'0 



(1 - K 2 f 



Q(k) 



3 2 tf - 



■ sin 2 § 



(26) 



with k = R 2 /ror c . In Eq. (|26l) Q(k) represents the mag- 
netic energy of the dipole field in a sphere having the ra- 
dius i?» = R 2 /r c < R and enclosed in the plasma sphere. 

Next, we evaluate the energy loss Wj(t) by the plasma 
which is determined by the surface current density, j local- 
ized within thin skin layer, R — 5<r<R + 5 with 6 — > 0, 
near plasma boundary. Therefore in Eq. (1221) the volume 
J? c is in fact replaced by the volume f2$ ~ R 2 8 of the shell 
with R — 5<r<R + 5. Within the skin layer we take 
into account that E = — -[v x H] and H r (R) = 0. Also 
it should be mentioned that the boundary of the volume 
f2$ moves with a velocity v and the electrical field has a 
jump across the plasma surface. The detailed calculations 
are given in [3] and we provide here only the final result 
which reads 



Qj(t) = - f j • Erfr = ^ / 

Jn s 7 Js R 



H 2 {R) 
8tt 



dS. 



(27) 



Here 7~ 2 = 1 — /3 2 and Sr are the relativistic factor and 
the surface of the expanding plasma, respectively. Note 
that the moving boundary modifies the surface current 
which is now proportional to the factor -f~ 2 [3"l[2"0"]. Equa- 
tion (|2"T|) shows that the energy loss by the plasma per unit 
time is equal to the work performed by the plasma on the 
external electromagnetic pressure. As shown in Ref. [1] 
this external pressure is formed by the difference between 
magnetic and electric pressures, i.e., the induced electric 
field tends to decrease the force acting on the expanding 
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Fig. 1. (Color online) The energy components Ws(t) (the lines without symbols) and W.j(t) (the lines with symbols) normalized 
to the quantity p 2 jr% for two values of the factor f3 as a function oict/ro calculated from expressions (|24|l and ()28[) with r c = 0.5ro- 
The solid and dashed lines correspond to the dipole orientations -& — and § = 7r/2, respectively. 



plasma surface. The total energy loss by the plasma sphere 
is calculated as 



Wj(t) = / Qj(t')dt 



2~4 (21 + i 



• I I — - (28) 



where £ = R/r$. In a non-relativistic case Eq. 
(see Appendix): 



Wj{t) 



e 



r* (l _ ?y 



cos 2 + 



1 + 



sin 2 i? 



yields 



(29) 



The change of the electromagnetic energy in a control 
sphere is calculated from Eq. (|2U)) . At R < r c < cf (the 
electromagnetic pulse fills the whole control sphere) it is 
straightforward to show that Z\Wem(*) = Ws(t) — Wj{t) 
as predicted by the energy balance equation ([2~Tj) . The 
non-relativistic limit of AW-EM_{t) can be evaluated from 
Eqs. ((26)) and ([29| . Also it should be emphasized that in 
Eqs. (pM|) . (|2"8)l as well as in Z\Wem(*) the dependence on 
the orientation i? of the dipole moment is determined by 
the factor ai (■&) given by (|2"5")) . At parallel (or antiparallel) 
orientation with $ = (or $ = 7r) this factor is ai — 
I + 1 and Eqs. (|2l)). (l28| as well as AWem^) coincide 
with the corresponding formulas derived in Ref. [1] . As in 
Eq. (0) it is convenient to represent the energy components 
W s (t) and Wj(t) in the form W s ,j(t) = W s ,j\\(t) cos 2 + 
Ws,j_i(t) sin 2 #, where W\\(t) and Wj_(t) are the energy 
components at parallel (d = 0) and transversal (d = n/2) 
orientations of the dipole moment, respectively. 

As an example in Figure[T]we show the results of model 
calculations for the time evolution of the energy compo- 
nents Wg (t) and Wj (t) for two extreme orientations of the 
dipole magnetic field with § = and d = tt/2. We recall 
that at ^ t ^ r c /c, i.e. the electromagnetic pulse does 



not yet reach to the surface of the control sphere, the en- 
ergy flux vanishes and Wg(t) = 0. Thus, for a given radius 
r c of the control sphere the energy flux occurs within the 
time interval At = (r c /c) (1//3— 1). Note that this interval 
decreases with the expansion velocity. As expected in Fig- 
ure [1] it is seen that both components of the energy Ws(t) 
and Wj (t) increase monotonically with time and strongly 
depend on the orientation d of the magnetic dipole. In 
addition, since the dipole magnetic field at the location 
of the plasma cloud is stronger at smaller orientation an- 
gles the energy flux Ws(t) and the energy loss Wj{t) 
are maximal for the parallel (or antiparallel) orientation 
of p. It is also noteworthy the dependence of the energy 
components on the relativistic factor f3 of the plasma ex- 
pansion. In the case of the non-relativistic expansion the 
energy flux Ws (t) is systematically larger than the energy 
loss Wj(t) (see FigureQ] left panel) at least in the wide in- 
terval of the expansion time (for comparison see Figure [2] 
below). Furthermore, both energies increase with the fac- 
tor /? but now Wj(t) exceeds Ws(t) in the wide interval 
of the expansion time (Figure [TJ right panel) . 

To gain more insight in Figure [2] we show also the 
time evolution of the ratios Is(t) = Ws(t)/Qo(t) and 
rj(t) = Wj(t)/Q (t) for ■& = and i9 = tt/2. The rela- 
tivistic factor f3 varies from the non-relativistic (f3 = 0.01) 
to the relativistic ((3 = 0.9) values. Here Qo(t) = Q(£) 
is the dipole magnetic energy escaped from the plasma 
sphere. Previously it has been found [15] that for the 
non-relativistic ((3 <C 1) expansion of a one-dimensional 
plasma slab and for uniform external magnetic field ap- 
proximately the half of the outgoing energy is gained from 
the plasma, while the other half is gained from the mag- 
netic energy, Ws — 2Wj ~ 2AWem- In the case of the 
expansion of a spherical plasma with (3 <C 1 in the uni- 
form magnetic field, Wj = 1.5Qo, ^W-em = 0.5Qo, and 
W s = 2Q , where Q = H$R 3 /6 [IS]. Therefore in this 
case the released electromagnetic energy is mainly gained 
from the plasma. 
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Fig. 2. (Color online) The ratios i"s(£) (the lines without symbols) and Pj(t) (the lines with symbols) for four values of the 
factor ft as a function of ct/ro calculated from expressions ()24|l and (|28[) with r c — 0.5ro. The solid and dashed lines correspond 
to the dipole orientations 1? = and i? = n/2, respectively. 



In contrast to the cases with uniform magnetic field 
discussed above in the present context there are no simple 
relations between the energy components Ws{t), Wj{t) 
and Qo(t). However, at the initial stage (t <C r c /v) of non- 
relativistic expansion the dipole field at large distances 
can be treated as uniform and the energies Ws(t) and 
Wj(t) are close to the values 2Q (t) and 1.5Qo(i) (Fig- 
ure [5]), respectively. It is noteworthy that for arbitrary 
[3 the ratios Fsit) and 2j(i) are only weakly sensitive 
to the orientation of the dipole magnetic field although 
the quantities Ws(t), Wj(t) and Qo(t) strongly depend 
on (see Figure Q] and Eq. © with (A.3), (A.4)). More- 
over, for any /3 the ratio -Tj(i) is almost constant and 
may be approximated as rj(i) ~ Tj(0) or alternatively 
Wj(t) ~ 1.5C*Q (i), where C* = 7 - 6 (l-^)- 4 (l + 2^)- 2 is 
some kinematic factor which is independent on the dipole 
orientation z9. These features merely indicate that for arbi- 
trary /3 the energy loss Wj(t) increases proportionally to 
the escaped energy Qo{t) while the parallel (Ws\\(t)) and 
the transversal (Wgj_(£)) components of the emitted en- 
ergy increase proportionally to the parallel (<3o||(*)) an d 
transversal {Qo±(t)) components of the escaped energy, 
respectively, with W S ||(t)/Q ||(<) a W s± {t)/Q 0± (t) ~ r s (-, 
Finally, in Figures[T]and[2]it is seen that for non-relativistic 



expansion the emitted energy is gained from both the 
plasma cloud and electromagnetic energy in the control 
sphere. However, at the final stage (t = r c /v) of the rel- 
ativistic expansion (with /3 ~ 1) Ws — Wj and hence in 
this case the emitted energy Ws is mainly gained from the 
plasma. 



4 Conclusion 

An exact solution of the uniform radial expansion of a neu- 
tral, infinitely conducting plasma sphere in the presence of 
a dipole magnetic field has been obtained for an arbitrary 
orientation of the dipole moment in the space. Thus we 
have generalized our result obtained in Ref. [J] for a paral- 
lel (or antiparallel) orientation of the dipole moment when 
p || ro. The electromagnetic fields are derived by using the 
appropriate initial and boundary conditions, Eqs. (|T5|) and 
(TT4"]) . respectively. As expected the electromagnetic fields 
are perturbed only within the domain located between 
the surfaces of the expanding plasma (r = R) and the 
■ information (r = ct) spheres. Outside the sphere r = ct 
the magnetic field is not perturbed and coincides with the 
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dipole magnetic field. In addition, since the dipole mag- 
netic field at the location of the plasma cloud is stronger 
for parallel (or antiparallel) orientation the induced elec- 
tromagnetic fields are maximal for p || ro. We have also 
studied the energy balance during the plasma sphere ex- 
pansion with an arbitrary orientation of the dipole. The 
model calculations demonstrate that the energy compo- 
nents increase monotonically with time and strongly de- 
pend on the orientation of the dipole moment (Figured]), 
being maximal at $ = 0,7r. On the other hand, the ratios 
rs(t) and -Tj(i) are only weakly sensitive to the variation 
of i9, the latter being a slowly varying function of time 
(Figure [2]). This is because the energy loss Wj(t) varies 
almost linearly with the escaped energy Qo(t) while the 
parallel (Wgi|(t)) and the transversal (Wsx(t)) compo- 
nents of the emitted energy vary linearly with the parallel 
(Qo\\(t)) and transversal (Qo±(t)) components of the es- 
caped energy, respectively. The calculations also show that 
for non-relativistic expansion the emitted energy is gained 
from both the plasma cloud and perturbations of the elec- 
tromagnetic energy in the control sphere. For relativistic 
expansion Ws — Wj and the emitted energy is practically 
gained only from the plasma sphere. 

Lastly, to solve the boundary value problem we have 
applied the conical flow technique which is only valid for 
the uniform expansion of the plasma cloud with constant 
v. In principle the solution for the arbitrary time-dependent 
velocity v(t) can be obtained by the Laplace transforma- 
tion technique used, for instance, in g]. This yields an 
integral equation for the vector potential which allows a 
general analytical solution only for a one-dimensional ex- 
pansion. For two- and three-dimensional expansions the 
general solutions are not known and these cases require 
separate investigations. 
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Appendix: Some summation formulas 

Using the known relation [3T] 



Similarly one obtains: 



(A.l) 



where x < 1, one can derive some summation formulas 
which are used in the main text of the paper. The first 
relation is evaluated as 



/i(*) = £ 



1(1 + 1)- 

2£+l 



„2f 



2x dx 



[*/(*)] 



8a; 2 - x 4 + 1 



1 . l + x 
m 



(1 - x 2 ) 3 2x ~l-x 

£ 2 (£ + l) 2e _ xd_ 
21 + 1 X ~ 2 dx 



(A.3) 



X 



x 4 + 8x 2 



1 



(1-x 2 ) 3 



1 . 1 + X 
— m 



2x 1 — x 



(A.4) 
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2\2 



1, 1 + X 
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In (A.2) we have used the boundary condition /(0) = 0. 



